Self-assembly of 2D membranes from mixtures of hard rods and depleting polymers 
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We elucidate the molecular forces leading to assembly of two dimensional membrane-like structures com- 
posed of a one rod-length thick monolayer of aligned rods from an immiscible suspension of hard rods and 
depleting polymers. We perform simulations which predict that monolayer membranes are thermodynamically 
stable above a critical rod aspect ratio and below a critical depletion interaction length scale. Outside of these 
conditions alternative structures such as stacked smectic columns or nematic droplets are thermodynamically 
stable. These predictions are confirmed using an experimental model system of virus rod-like molecules and 
non-adsorbing polymer. Our work demonstrates that collective molecular protrusion fluctuations alone are suf- 
ficient to stabilize membranes composed of homogenous rods with simple excluded volume interactions. 



Colloidal membranes are two dimensional (2D) surfaces 
composed of a one rod-length thick monolayer of aligned 
nanorods. Equilibrium formation of such structures requires 
assembly to readily propagate in two dimensions while self- 
limiting the third. Previous approaches towards assembly 
of colloidal membranes utilized chemically heterogeneous 
rods that mimic the dichotomic structures of amphiphilic 
molecules comprising biological membranes [1]. Here, we 
use computer simulations to demonstrate that structurally 
and chemically homogeneous hard rods can form equilibrium 
monolayers in the presence of depletant molecules, suggest- 
ing that geometry as well as chemical heterogeneity can be 
used to design assembly pathways of self-limited structures. 
Furthermore, we discover bounds on the molecular parame- 
ters that support formation of equilibrium membranes. These 
results have fundamental as well as practical significance. Ex- 
tensive research has shown that hard particle fluids undergo 
entropy-driven assembly into a myriad of 3D structures IeHSD- 
Our work demonstrates that entropic forces can also drive for- 
mation of 2D structures. From a practical perspective, equi- 
librium colloidal membranes may enable manufacture of in- 
expensive and easily scalable optoelectronic devices |@]. 

Our study is motivated by recent experiments on suspen- 
sions of monodisperse rod-like colloidal viruses and the non- 
adsorbing polymer Dextran [7] (Fig. [T). fd viruses alone 
approximate the behavior of homogenous rods interacting 
with repulsive hard-core interactions [8]. The polymer in- 
duces an entropy-driven attractive (depletion) potential be- 
tween the rods, the strength and range of which can be tuned 
by changing the polymer concentration and radius of gyra- 
tion respectively (Fig UK) I2D. At high polymer concentra- 
tions viruses condense into smectic-like stacks of 2D mem- 
branes (FigHt>) m. With decreasing polymer concentra- 
tion (attraction strength), individual 2D monolayers (mem- 
branes) within a smectic filament unbind, indicating that the 
membrane-membrane interaction switches from attractive to 
repulsive J7|] (FigHt). 

Understanding the molecular origin of the repulsive 
membrane-membrane interactions is our primary goal. Ex- 
periments revealed significant protrusions of rods from iso- 
lated colloidal membranes, the magnitude of which could be 
tuned by changing the concentration of non-adsorbing poly- 



mer |7|]. In contrast, these fluctuations were suppressed in 
stacked membranes. It was proposed that the entropy penalty 
associated with suppressing protrusion fluctuations of individ- 
ual rods leads to repulsive interactions that stabilize isolated 
membranes under moderate osmotic pressure [11]. However, 
other plausible factors could also stabilize membranes, includ- 
ing attractive interactions between virus tips and depletant, re- 
pulsions due to bending (Helfrich) modes, or kinetic trapping 
of membrane intermediates. To elucidate these issues, we de- 
velop a computational model which demonstrates that protru- 
sion interactions alone are sufficient to stabilize membranes in 
equilibrium. In contrast to the previous model which consid- 
ered only protrusions of isolated rods [7], our work indicates 
that collective protrusion undulations dominate repulsive in- 
teractions between membranes. Surprisingly, the simulations 
predict that membranes are stable only for a certain range of 
rod aspect ratios and depletant sizes. We experimentally con- 
firm the latter prediction. 

We model the fd rods as hard spherocy finders with diam- 
eter a and length L. The non-adsorbing polymer is mod- 
eled as ghost spheres of diameter 6 which freely interpene- 
trate one another but behave as hard spheres when interacting 
with rods [12]. Compared with an effective pair potential ap- 
proach, this model accounts for multi-rod interactions induced 
by polymers Jl3r - [l5ll . We perform Metropolis Monte Carlo 
(MC) with periodic boundary conditions fll6l1 . The total num- 
ber of rods 7V r is fixed, the sphere osmotic pressure p s is set by 
insertion/deletion moves, and constant pressure is maintained 
in the xy plane by performing volume-change moves, while 
the box size is fixed in z direction 111711 . Simulation results 
are reported with a as the unit of length, feT as the unit of 
energy, and ksTa~ 3 as the unit of pressure. 

Membrane-like structures have two generic repulsive in- 
teractions of distinct origin which dominate at different sep- 
aration lengthscales. At large separations, slowly decaying 
low energy bending (Helfrich) modes dominate 111 811 . In con- 
trast, at separations comparable to the monolayer thickness 
(the rod length), protrusions of molecules from the membrane 
surface generate the primary repulsive force [H9I . We expect 
that long ranged bending modes provide a negligible contri- 
bution to the stabilization of colloidal membranes for the fol- 
lowing reasons. First, the range of the depletion attraction 
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FIG. 1. Schematic illustrations and optical micrographs of the self- 
assembled structures observed in suspensions of the filamentous 
virus fd and non-adsorbing polymer |7]. A) Non-adsorbing poly- 
mer induces effective attractive interactions between rods. B) DIC 
micrograph and schematic of a nematic tactoid formed at low deple- 
tant concentration. C) At intermediate depletant concentrations, rod- 
like viruses condense into macroscopic one rod-length 2D fluid-like 
membranes. D) At high depeltant concentration, membranes stack 
on top of one another, forming smectic filaments. All scale bars are 
5/xm. 



that balances repulsive forces to drive membrane stacking is 
comparable to the depletant size. On this scale protrusion re- 
pulsions dominate. Second, the bending modes involve de- 
viations of rods from the preferred direction and their inter- 
action strength scales with bending modulus as /bend ~ Kq 1 
lEoi Elll . The large bending modulus measured for fd mem- 
branes k c = 150/cbT thus results in very weak Helfrich repul- 
sions Hill]]. Consistent with these arguments, a theoretical 
calculation (SI Fig. 5) shows that the protusion interactions in 
simulated membranes exceed the strength of the Helfrich in- 
teractions by four orders of magnitude at relevant separations. 

Based on the preeminence of protrusion modes, which do 
not involve rod tilting, in most simulations we restrict sphe- 
rocylinder orientations to be perfectly aligned along the z di- 
rection. This simplification greatly enhances computational 
efficiency, allowing us to extensively map the phase diagram 
as a function of all relevant molecular parameters. Our ap- 
proximation is justified by Fig. 5 in the SI and the fact that 
simulations in which the fixed orientation constraint is relaxed 
predict similar phase behavior and membrane-membrane in- 
teractions (e.g. Fig.0. 

Membrane-membrane interaction potential We first use 
umbrella sampling [16] to measure the free energy per rod 
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FIG. 2. (left) Free energy per rod of an interacting membrane pair, 
f(d) — /(oo), plotted as a function of membrane surface separation 
d — L, shown for three different depletant concentrations indicated 
by values of p s , with sphere diameter S = 1.5 and aspect ratio L = 
100. The dashed line is the free energy calculated with orientational 
fluctuations atp s = 0.06. (right) Snapshots of two membranes from 
unbiased trajectories, (top) Membranes attract at p s — 0.12. (middle, 
bottom) Snapshots for p s = 0.06 from simulations with (middle) 
parallel rods and (bottom) rods with orientation fluctuations. In both 
cases membranes drift apart, indicative of a repulsive potential. 

/, as a function of the separation between the centers of mass 
of two membranes, d (Fig.0. At low osmotic pressures (e.g. 
p s = 0.06), f(d) — /(oo) has no attractive region sufficient 
to overcome translational entropy; i.e., the stacking of disks is 
suppressed and the isolated colloidal membrane phase is sta- 
ble. For larger osmotic pressures (p s > 0.08), the free energy 
has a substantial minimum at finite membrane separations, 
signifying that membranes will stack to form the smectic-like 
columns. Consistent with these free energy results, unbiased 
simulations for these parameters resulted in two membranes 
which were respectively isolated and stacked at low and high 
osmotic pressures, as shown in Fig. [2] (right). The free en- 
ergy and a representative snapshot are also shown for rods 
with orientational fluctuations at p s = 0.06. Note that isolated 
membranes are stable and the interaction free energy is com- 
parable to the case with parallel rods; the repulsion is slightly 
weaker with orientational fluctuations because they decrease 
the equilibrium areal rod density. 

Phase diagram. We computed the equilibrium phase be- 
havior as a function of osmotic pressure, rod aspect ratio, and 
sphere diameter as follows (Fig. [3}. To identify the nematic - 
membrane phase boundary, we performed separate unbiased 
simulations starting from initial conditions in which (1) rods 
have random positions and (2) rods are aligned in a flat layer. 
For all results shown, the simulation outcomes were inde- 
pendent of initial conditions. To identify the transition from 
membranes to smectic filaments, a parameter set was consid- 
ered to yield smectic layers if the total free energy of the 
attractive basin in the membrane-membrane interaction po- 
tential satisfies F < Fo = k^,T\iip m VQ with exp(— /3F) = 
ff(s)<o dsex.p(— 2(3Mf(s)), M the number of rods in one 
membrane, vq a standard state volume, and p m a membrane 
concentration. A finite value of Fo accounts for membrane 
translational entropy. We roughly estimate M = 10 4 and 
PmVo = 10 -8 from the experimental conditions; the location 
of the phase boundary is not sensitive to the value of p m vq . 

Fig. [3]A illustrates the location of the equilibrium nematic 
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FIG. 3. Phase diagrams from simulation and experiment. Triangles A denote denote parameters that lead to nematic configurations, + symbols 
correspond to isolated membranes, and ■ symbols correspond to smectic layers. (A),(B) Phase diagrams determined from simulations for 
varying osmotic pressure p s and (A) aspect ratio L with sphere diameter (polymer radius of gyration) S = 1.5, and (B) varying sphere 
diameter with L = 100. The solid lines identify the the isolated membrane/smectic and nematic/isolated membrane phase boundaries. They 
are fit by eye to simulation results except for the nematic/isolated membrane boundary in (A), which is a theoretical prediction Il22ll . (C) The 
experimental phase diagram corresponding to (B) using mixtures of fd viruses and PEG/PEO polymers. The final concentration of viruses 
was fixed at 5mg/mL and both polymer concentration and molecular weight were varied to change osmotic pressure and polymer radius of 
gyration, R g , respectively. As noted in the text, chiral structures such as helical ribbons which appear near the nematic and isolated membrane 
boundary are not shown. 



phase, isolated membranes, and smectic stacks as a function 
of rod aspect ratio and depletant concentration. Interestingly, 
isolated membranes are thermodynamically stable over a sig- 
nificant span of osmotic pressures, but only for rods with as- 
pect ratios larger than L = 30. Simulations with orientational 
fluctuations also indicate a minimum aspect ratio for stable 
membranes, which is somewhat larger. These predictions are 
consistent with previous simulations of rods with L = 5 that 
did not find equilibrium monolayers fll4l1 . The disappearance 
of the isolated membrane phase for shorter rods arises from 
the interplay between the geometry of rod-like particles and 
attractive depletion interactions. Since the strength of the at- 
tractive interaction between two rods scales linearly with rod 
length, increasing the rod length lowers the osmotic pressure 
associated with the nematic to membrane transition. On the 
other hand, the transition from isolated membranes to smectic 
filaments is determined by the roughness of colloidal mem- 
branes, which is independent of rod length but decreases with 
increasing depletant concentration. Based on this argument, 
the location of the transition between colloidal membranes 
and smectic filaments should be independent of rod length, 
which is indeed observed for rod lengths between 30 and 
100. For longer rods the location of the transition slightly de- 
creases with increasing rod length, due to 2D crystallization 
of rods within membranes (see the SI for details of membrane 
crystallization and a determination that finite size effects do 
not affect the results). At a critical rod length the nematic - 
membrane phase boundary intersects the membrane- smectic 
filament phase boundary, ending the equilibrium membrane 
phase. 

Fig. [3j3 reveals that the depletant size 6 also influences the 
topology of the phase diagram. For S > 1.7 colloidal mem- 
branes are unstable at all osmotic pressures and there is a di- 
rect transition from the nematic phase to smectic filaments. 



In contrast, for 5 < 1.7 colloidal membranes are the equilib- 
rium phase at intermediate depletant concentrations between 
a low osmotic pressure nematic phase and high osmotic pres- 
sure smectic filament phase. Decreasing the depletant size fur- 
ther below this critical value significantly expands the range 
of osmotic pressures for which colloidal membranes are sta- 
ble. These results can be understood as follows. Increasing 
the depletant size expands the effective range of the attractive 
potential between two colloidal membranes, which in turn re- 
quires longer range repulsive interactions to stabilize colloidal 
membranes. For large enough depletant molecules, the repul- 
sive protrusion interactions are not sufficiently long-ranged to 
overcome the attractive potential and colloidal membranes be- 
come unstable for all osmotic pressures. 

Origins of monolayer stability. To understand the nature 
of the repulsive membrane-membrane interactions, we de- 
termine their functional form / pr by subtracting the deple- 
tion interaction fa from the measured membrane-membrane 
free energy, f pT (d s ) = f(d) — fd(d). The depletion term 
is given by fd(d) = Ps(v QX }d s , where v QX is the volume 
excluded to spheres by rods, and (-)d indicates an ensem- 
ble average over configurations at a particular separation 
d. We then adapt a calculation in Ref. [21] to obtain the 
membrane-membrane interaction due to collective protrusions 
as 2^2d/pr = exp[— irj(d — L) 2 /3/cbT] with 7 the surface 
tension, p2d the area per rod, and B a constant. The calcu- 
lations are presented in further detail in the SI. As shown in 
Fig. |4l the measured repulsive interaction / pr is well described 
by this functional form, with fit values of 7 that are close to the 
surface tension extracted from simulated height-height corre- 
lation spectra (SI Fig. 1). Thus, the membrane-membrane 
repulsion primarily arises from collective protrusion undula- 
tions. 

Experimental phase diagram. Simulations predict a crit- 
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FIG. 4. The protrusion interaction potential is well-fit by the theory 
in some parameter ranges. The dotted lines show the repulsive in- 
teraction potential f pr measured from simulations and the solid lines 
correspond to the best fit to the protrusion undulation potential given 
in the text with B and 7 as fit parameters. Parameters are L = 100, 
S = 1.5 and (A) p s = 0.06, (B) p s = 0.08 and the best fit values are 
(A) B = 0.8, 7- 1 = 213, (B) B = 0.9, 7 _1 = 156. 



tant size. While most previous simulations of hard rods con- 
sidered small aspect ratios, our prediction of a critical aspect 
ratio below which the colloidal membrane phase disappears 
suggests that large aspect ratios are crucial for the phase be- 
havior observed in Ref . Q] . The predicted critical aspect ratio 
is only qualitative, but can be tested by monitoring the phase 
behavior of depletant and rods with varying lengths, as the 
prediction of a critical depletant size was tested here. 
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ical depletant size above which isolated membranes are un- 
stable with respect to stacks of membranes for all osmotic 
pressures. We experimentally verify this prediction using a 
mixture of fd virus and non-adsorbing polymers(for methods 
see the SI). As shown in [3j there is qualitative agreement 
between simulations and experiments in two respects. First, 
colloidal membranes are unstable for depleting polymer of 
large size; i.e. there is a direct transition from the nematic 
phase to smectic filaments. In contrast, for smaller poly- 
mer sizes, colloidal membranes are stable. Second, with de- 
creasing polymer size the osmotic pressure (polymer concen- 
tration) at the transition from colloidal membranes to smec- 
tic filaments increases. Several points need to be considered 
when comparing the experimental and computational phase 
diagrams. First, there is a gap in the data between the poly- 
mer sizes corresponding to R g = 9.7nm and R g = 17.9nra 
due to limited commercial availability of polymers with ap- 
propriate size. Second, the transition pressure from the ne- 
matic/isotropic phase to colloidal membranes increases pre- 
cipitously for smaller polymer sizes (R g < 5.2 nm). This is 
due to the deviations of the/d system from an ideal model hard 
rod system due to its surface charge. Making the depleting 
polymer size smaller than the electrostatic repulsion length 
greatly reduces the strength of the attractive interactions, re- 
quiring a higher depletant concentration to induce condensa- 
tion of colloidal membranes J23I1 . Third, while the chirality of 
the individual viruses can influence the assembly pathways, 
we have determined that the locations of transitions in the ex- 
perimental phase diagram are independent of the chirality of 
the constituent rods. 

In summary, this study demonstrates for the first time that 
entropic forces are sufficient to stabilize monolayer colloidal 
membranes at equilibrium. We find that collective protrusion 
undulations are the primary force that stabilizes isolated mem- 
branes. While experimental observations of protrusions in 
Ref. 10] were inferred as individual rods, those experiments 
only fluorescently labeled a small fraction of rods and thus 
could not resolve collective modes. 

The simulations also predict that the width of the isolated 
membrane phase depends strongly on aspect ratio and deple- 
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